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We investigate the spin fluctuations and the pairing symmetry in A :E Fe2-j,Se2 by the fluctuation 
exchange approximation. Besides the on-site interactions, the next-nearest-neighbor antiferromag- 
netic coupling J2 is also included. We find that both the itinerant and the localized natures of 
electrons are important to describe the recent experimental results of the spin fluctuations and 
the pairing symmetry. In particular, a small J2 coupling can change the pairing gap from the d- 
wave symmetry to the s-wave symmetry. We have also studied the real-space structures of the gap 
functions for different orbits in order to gain more insight on the nature of the pairing mechanism. 
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I. INTRODUCTION 

The discovery of a new family of iron-based supercon- 
ductors (FeSC) A x Fe 2 - y Se 2 (A=K, Rb, Cs) with T c more 
than 30 K has attracted much attention 1 - - — . The experi- 
ments on these materials show many surprising results 
which apparently differentiate them from other FeSC 
families (iron chalcogenide and iron pnictide). First, 
angle-resolved photoemission (ARPES) studies indicate 
that there are only electron-like Fermi surfaces (FS) and 
no hole-like FS in A a .Fe 2 _j / Se^— , while other FeSC fami- 
lies have both hole- and electron-like FS simultaneously 8 . 
Second, they are the first FeSC that approximates an 
antiferromagnetic (AFM) insulator accompanied by Fe- 
vacancy order 4 . And, the AFM order has a novel -\/5 x y/E 
block-type structure with an unprecedented high transi- 
tion temperature of Tjy = 559K and a large magnetic mo- 
ment of 3.31/zs/Fe^. The implications of these anoma- 
lous features on superconductivity in A x Fe2- a Se2 are not 
clear presently. And moreover, the questions about the 
pairing mechanism and the gap structure of the super- 
conductivity in these systems still remain open. 

From the theoretical viewpoint, the absence of hole 
Fermi pockets formally violates the nesting condition of 
the s± pairing symmetry, whose gap function changes 
signs between the hole and the electron pockets as in the 
most cases of other FeSC families, so that it requires a re- 
consideration of the itinerant spin- fluctuation mechanism 
in these new systems. Based on the multi-orbital Hub- 
bard model and the functional-renormalization-group 
calculation, Wang et alr^ find that the d x 2_ y 2-w&ve in- 
stability is the leading pairing channel with a subleading 
extended s-wave component. A similar conclusion has 
also been obtained by Maier et alM through performing 
the random-phase-approximation (RPA) calculation. On 
the other hand, based on the local antiferromagnetic ex- 
change interactions, Seo et al^ predict an s-wave pairing 
symmetry in this system that can account for the ARPES 
experimental results^—. Khodas et al. 13 suggest an s" 1 
pairing symmetry, in which the gap changes sign between 
the hybridized pockets, by including the interpocket pair- 



ing. 

Experimentally, the ARPES measurements have re- 
ported a nodeless superconducting (SC) gap on the 
large Fermi pockets around the zone corner in these 
mat er ials^~— ' 1 4 1 1 5 . In particular, the recent measurements 
further find an isotropic SC gap distribution on the small 
electron Fermi pocket around the Z poin t 16 ' 17 , which fa- 
vors the s-wave pairing symmetry. The experimental ev- 
idences from the nuclear magnetic resonance (NMR) ex- 
periments suggest an s" 1 -wave pairing symmetr y 18 ' 19 , 
which is consistent with the isotropic gaps reported by 
ARPES. And, the muon spin spectroscopy experimentSS 
suggest an s-wave symmetry and the results are incom- 
patible with a clean d-wave model. In addition, the 
specific heat measurement also indicates that the gap 
is nodeless^i. Though, the inelastic neutron scattering 
(INS) measurement o 22 ' 23 report a spin resonance in the 
superconducting state at Q w (7r,7r/2) which is very 
close to the value predicted by Maier et alr^- with a 
d-wave pairing symmetry, another spin excitation near 
Q = (0, 7r) was also observed in the superconducting 
state 23 . The latter is identical to that in many iron 
arsenide superconductors, which are believed to be the 
-wave superconductors. 

In this paper, we will study the spin fluctuations and 
the paring symmetry in A.;;Fe2-j,Se2 with the fluctuation 
exchange (FLEX) approximation by including the next- 
nearest-neighbor (NNN) antiferromagnetic (AF) cou- 
pling J2 , a localized magnetic interaction, into the multi- 
orbital Hubbard model which is treated within an itin- 
erant electron picture. In fact, the INS experiments 2 - 4 ^— 
have indicated that the NNN exchange coupling is AF 
and has a similar magnitude in many FeSC, although the 
nearest-neighbor (NN) exchange coupling is quite differ- 
ent. This implies that the NNN exchange coupling may 
play an important role in the FeSC. But, in the FeSC with 
both hole- and electron-like FS the NNN exchange cou- 
pling has the similar effect with the FS nesting on the spin 
excitations. However, we will show that both the itiner- 
ant and the localized natures of electrons are essential to 
describe properly the physical properties in A^Fe^-ySe^, 
where the hole-like Fermi pocket disappears. In partic- 
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ular, a small J2 can change the pairing gap from the 
d-wave symmetry to the s-wave symmetry, which is con- 
sistent with nodeless gap observed in recent experiments. 
Considering the degeneracy of the two pairing states near 
the transition, we suggest a possible s + id-wave state in 
order to account for the spin resonance observed in INS 
experiments. We have also studied the real-space struc- 
tures of the pairing gap at different orbits, from which 
we obtain the expressions for the favorable pairing gap 
functions. 



II. MODEL AND FLEX METHOD 

The model Hamiltonian consists of two parts: 

H = H + H int . (1) 

The non-interacting part Hq is described by the two- 
dimensional five-orbital tight-binding model in the un- 
folded (1-Fe) Brillouin zone (BZ) as introduced in Reffiol. 
The parameters of the model involves the nearest- and 
second-neighbor hopping among the Fe 3d orbits. The 
resulting tight-binding Hamiltonian reads 

5 

#o=]T ]T c{ a(y K aP (k)c kPa , (2) 

k,a a,/3— 1 

where c* kaa (cfc Qcr ) creates (annihilates) an electron with 
spin a and momentum k in the orbital a. The parame- 
ters used in constructing K a p(k) can be found in ReflTol. 
However, in this paper, we use a gauge transformation 
on the five orbitals to change the orbital symmetry from 
(dz 2 ,dxz,dYZ,dx 2 -Y 2 ,dxY) to (d z 2 : d xz ,d yz ,d xy ,d x 2_ y 2), 
where x, y refer to the nearest-neighbor Fe-Fe directions 
and X, Y refer to the diagonal directions. The transfor- 
mation matrix is given by 

/ 1 \ 

75 "75 

G = 75 75 ■ ( 3 ) 
1 
V -1 / 

and K a fi(k) can be obtained from H a p{k) in Reffiol 
through the following transformation: 

K{k) = G^H{k)G. (4) 

The band structure and the associated FS for 0.3 elec- 
tron doping (6 electrons per Fe for parent compound) 
are shown in Fig. [TJ We find that there are only electron 
FS around (0,±tt) and (±7r, 0), which is consistent with 
the results of ARPES experiments^—. Here, we do not 
consider the influence of Fe vacancies, because many ex- 
periments suggest that the antiferromagnetism with Fe 
vacancies and the superconductivity without Fe vacan- 
cies are in fact spatially separated^— . 




FIG. 1. (Color online) Energy band structure of the five or- 
bital model (a) and the FS in the unfolded BZ for 0.3 electron 
doping (b). The colors indicate the majority orbital character 
(red: d xz , green: d yz , blue: d xy , yellow: d x 2_ y 2, pink: d z 2, 
orange: d xz and d yz with the same weight). The arrows with 
solid lines in (b) indicate the weak nesting vectors q = (n, n) 
and (— 7r, — 7r), the arrows with dashed lines indicate the weak 
nesting vectors q = (71-, 0.67r) and (— 0.67r,7r). 

The interactions between electrons are included in Hi nt 
as following, 

Hint =U*^2 n 'iat n ial + V ^ n iarT THfi a ' 

— Jh ^ Si a ■ Sip + J2 ^ Si a ■ Sj/3, (5) 

where U (V) is the intra-orbital(inter-orbital) Coulomb 
interaction, Jh the Hund's coupling, and J2 the NNN 
antiferromagnetic coupling. n,i aa — c\ aa Ci acr and Si a ~ 

Yla-a-' c \aa T °°' c ia°' are * ne density and spin operators re- 
spectively, where r are the Pauli matrices. 

As mentioned above, due to the absence of the hole 
pockets in A^Fe2- y Se2, the NNN antiferromagnetic cou- 
pling J2 may play an essential role in determining the 
physical properties in these materials. Moreover, since 
the J2 exchange coupling originates mostly from superex- 
change processes through Se, the effect of J2 is under- 
estimated in the itinerant-electron model based on the 
pure iron lattice with only the on-site interactions. On 
the other hand, the INS experiments^— have indicated 
that the NNN exchange couplings are AF and have the 
similar values in many FeSC. Hence, we consider the ef- 
fect of an admixture of both the on-site interactions and 
the J2 exchange coupling here. A study based on the sim- 
ilar consideration for the interactions in the two-orbital 
model for pnictide superconductors with both electron 
and hole FS has been carried out from the strong cou- 
pling expansion^. 

We carry out the investigation using the FLEX 
approximation^— , in which the Green's function 
and the spin/charge fluctuations are determined self- 
consistently. Here, we consider the following spin fluc- 
tuations X s (<7, iuj n ) and charge fluctuations x c (q,iw n ): 

X s a p(q,i"n) = [ T dTe^(T T [S a (k,r) ■ Sp(-k,0)])> 
Jo 

(6) 
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% tfS (q,iu n ) = / dTe iul ™ T (T T [n a (k,T)np(—k, 0)]), (7) 
Jo 

where S(k, r) = e TH S(k)e~ rH with S(k) the Fourier 
transformation of Si, T is the temperature, r and w„ are 
respectively the imaginary time and the Matsubara fre- 
quency. For the five-orbital model, the Green's function 
G and the self-energy £ are expressed in a 5 x 5-matrix 
form. As we use the spin and charge fluctuations ex- 
pressed in Eq. © and Eq. ([7]), the susceptibility x° 
and the effective interaction V also have a 5 x 5-matrix 
form. The Green's function satisfies the Dyson equation 
G(fc) -1 = G°(fc) _1 — where the self-energy is given 

by S m „(fc) = lEgE^^.mi/ls)^!/^ - q) and the 

bare Green's function reads G°(k) — [iuj n — K(k) + /u] -1 . 
Here, k includes both momentum and frequency with 
k = (k,iuj n ). The fluctuation exchange interaction is 
given by: 

v(q) = §&(«)[(/ - x ^*?))- 1 - ip s (q) 
+ \Hq)[{i + x\q)UM)- 1 ~i]uM 
+ l(u s (q)x°(q)Us(q) + u c (q)x°(q)u c (q)l (8) 

with the identity matrix I and the irreducible susceptibil- 
ity X%(q) = ~|f Efc G 0a (k + q)G a p(k). The interaction 
matrix for the spin(charge) fluctuation U S (U C ) is given 
by: For a = /3, t/^, = U - 4 J 2 cosfc^ cosfc. y (U° p = U)\ 
For a ^ 0, U s a[j =J H - 4J 2 cos k x cos fc a (U c afj = 4V). 

The Dyson equation, the self-energy and the interac- 
tion matrix Eq. © form a closed set of equations and 
can be solved self-consistently. Since the equations must 
be solved at fixed electronic density n, they are subject 
to the additional constraint: 

n (T) = ^J2 ^"° + TrG(fc, fJ,(T)) = n, (9) 
k 

where the chemical potential fj, must be adjusted accord- 
ingly during the self-consistent calculation. 

After obtaining the renormalized Green's function G, 
we can solve the "Eliashberg" equation, 

\<f> af > « = - § E E V S ( fc - 1) 

q liv 

xGp u (q- k^uik- q), (fO) 

where the spin-singlet and spin-triplet pairing interac- 
tions V s and V t are given by, 

v s (q) = lu s r(q)u s -lu c r(q)u c +l(u c +u s ), (ii) 



v\q) = ~^ s x s (q)u s - ^u c x c (q)u c + ^(u c - u s )- 

(12) 



The most favorable SC pairing symmetry corresponds to 
the eigenvector 4> a p(k) with the largest eigenvalue A. 

Our numerical calculations are carried out on the 
64 x 64 k meshes with 4096 Matsubara frequencies. In 
the calculations, we set the intra-orbital Coulomb inter- 
action U — 1.2eV, the inter-orbital Coulomb interaction 
V = 0.35eV and the Hund's coupling J H = 0.2eV. The 
electron density is set as 6.3 electrons per Fe and the 
temperature is 0.008eV. 



III. RESULTS AND DISCUSSION 
A. Spin fluctuation 

As the spin fluctuations dominate over the charge fluc- 
tuations in our model described above, we only discuss 
the spin fluctuations and their contributions to the pair- 
ing interactions of electrons in this paper. We first 
study the evolution of the spin fluctuations with the 
NNN antiferromagnetic coupling J 2 . In Fig. [2J we 
present the physical spin susceptibilities x s (<Z,w = 0) = 
^ a(3 X^ /3 (q,w = 0) for four different values of J 2 . For 
,/ 2 = 0, the absence of the hole pocket around T point 
in the unfolded BZ reomves the dominant q = (0, ±7r) 
and (±7r, 0) nestings in iron pnictide, so that the peaks 
of the spin fluctuations at (0, ±tt) and (±tt, 0) disap- 
pear in A a: Fe 2 _j ; Se 2 . Instead, there are weak nestings at 
q = (±7r, ±0.67r), (±0.67r, ±71") and (±7r, ±7r) in this sys- 
tem (see Fig. DIb)), which lead to the broad peak around 
(it, 7r), together with the four satellite peaks at (tt, ±0.67r) 
and (±0.67r,7r) (Fig. Ufa)). We note that these results 
are not consistent with the experimental results from INS 
measurement o 22 i 23 , in which the peak around (tt, tt) are 
not observed. Thus, a simple Fermi surface nesting sce- 
nario within the itinerant-electron model is not enough to 
explain the neutron scattering experiments. When turn- 
ing on the J 2 interaction term, we find that the broad 
peak around (n, 7r) decreases rapidly. Meanwhile, the 
intensity of the four satellite peaks increases and their 
positions move towards (0,±7r) and (±7r, 0), as shown 
in Fig. [U (b)-(d) for J 2 = O.OleV, 0.04eV and 0.05eV 
respectively. As the J 2 term favors the diagonal antifer- 
romagnetic order in the square lattice of Fe, it competes 
with the spin fluctuations coming from the nesting of FS 
with the nesting wavevector (tt, tt) which is a kind of axial 
antiferromagnetic fluctuations. Thus, an enough large J 2 
can suppress the spin fluctuations from the nesting effect 
of the FS. It is worthy to pointing out that we do not need 
a large value of J 2 comparing to the band width (4eV) 
to change the character of the spin fluctuations. We find 
that when J 2 > 0.04eV, the peak at (tt, tt) is completely 
suppressed which is consistent with the results of INS 
experiment o 22 ! 23 . And, the (7r, 0) spin excitations ob- 
served in Rbo.82Fei.68Sesr^ also emerge for J 2 > 0.04eV. 
In addition, for 0.03eV < J 2 < 0.05eV, the peaks around 
(±7r, ±0.67r) and (±0.67r, ±7r) move to about (±7r, ±0.57r) 
and (±0.57r, ±7r), which is also consistent with the INS 
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FIG. 2. (Color online) Physical spin susceptibilities for dif- 
ferent NNN antiferromagnetic coupling J 2 . (a) J 2 = 0, (b) 
J 2 = O.OleV, (c) J 2 = 0.04eV and (d) J 2 = 0.05eV. 

results. If we further increase J2 to J2 > 0.05eV, the 
peaks around (±7t, ±0.57t) and (±0.57r, ±7t) which comes 
mainly from the FS nesting disappear completely, and 
the four peaks of the spin fluctuations will now at the 
positions between (jr, 0.57r) and (n, 0), and its symmetric 
points, so that it breaks down the consistence with the 
INS experiments. These results imply that an admixture 
of both the on-site interactions and the J2 exchange cou- 
pling is appropriate and the J2 exchange coupling has sig- 
nificant effects on the physical properties of A a; Fe2_7 / Se2. 
So, we will further study the pairing symmetries medi- 
ated by the spin fluctuations with the including of the J2 
term. 



B. Symmetry of the superconducting pairing 

With the spin fluctuations, we can search for the most 
favorable pairing symmetry mediated by the spin fluctu- 
ations by solving the Eliashberg equation Eq.(|10p. We 
find that the eigenvalue A has the maximum values in 
the spin-singlet channel and it is nearly zero in the spin- 
triplet channel in the range of parameters we have used. 
The paring functions corresponding to the two leading 
eigenvalues in the spin-singlet channel have the s-wave 
and d- wave symmetry. In Fig. [3J we show these eigen- 
values as a function of J 2 for 0.3 electron doping corre- 
sponding to K .6Fe2Se2. For a small J2 (< 0.045eV), the 
eigenvalue for the d wave is larger than that for the s- 
wave, so the ci-wave pairing symmetry is dominant over 
the s-wave pairing symmetry. However, the eigenvalue 
for the s wave pairing follows a monotonic increase with 
J2, and when J2 > 0.045eV the s-wave pairing surpasses 
the d-wave pairing and becomes the most favorable pair- 
ing state. We note that the threshold value for J2 is much 
smaller than that needed in the local-antiferromagnetic- 




0.00 0.02 0.04 0.06 
J 2 /eV 

FIG. 3. (Color online) J 2 dependence of the eigenvalues 
for the s-wave and the d-wave solutions of the Eliashberg 
equation at U = 1.2eV, V = 0.35eV, J H = 0.2eV and 
T = 0.008eV. 



exchange model 1 - and is comparable to the experimental 
value from the INS measurement) 2 ^. 

The gap functions on the FS for the s-wave and the d- 
wave symmetry corresponding to two representative val- 
ues of J2 = and 0.05eV are shown in Fig. 2) Here, 
the magnitudes of the gap have been normalized into the 
range from —1 to +1. In the case of J2 = 0, the s-wave 
gap function has nodes (zero points of gap) on the FS, 
yet the dominant pairing symmetry here is the d-wave. 
The gap of the d wave symmetry for J2 — [Fig. 0] 
(b)] has no node on the FS, this is consistent with the 
ARPES results of a nodeless SC gap on the large Fermi 
pockets around the zone corner—"— ' 14 ' 15 . But, it has line 
nodes along the diagonal direction in the Brillouin zone. 
Recent ARPES experiments further report an isotropic 
SC gap on the small electron Fermi pocket around the Z 
poin t 16 ' 17 , so it rules out the d-wave symmetry. In ad- 
dition, as argued by Mazin 38 , this d-wave gap is rather 
fragile when the hybridization between the electron pock- 
ets in the folded BZ is included. For the same reason, the 
d-w&ve gap for J2 = 0.05 [Fig. 2] (d)] is also fragile. The 
only pairing symmetry that is stable and consistent with 
experiments is the s-wave gap in the case of J2 = 0.05 
[Fig. 2] (c)], which is also the dominant pairing symme- 
try for J2 — 0.05 as shown in Fig. [3] We note that the 
gap magnitude of this s-wave state has in fact a varia- 
tion along the FS, though it has no node. To compare 
with the ARPES experiments, in Fig. 2] (e), we present 
the variation of the gap along the FS in the folded (2-Fe) 
BZ. We find that the gap on the inner FS is larger than 
that on the outer FS and it has a relative small variation 
(about 25%) along the FS. 

Though the nodeless s-wave gap is consistent with the 
ARPES experiments, it can not account for the spin 
resonance observed in the neutron scatterin g 22 ' 23 , be- 
cause a sign change between the gaps connecting by the 
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FIG. 4. (Color online) Pairing functions on FS. (a) and (b) are 
the s-wave and the d-wave solutions for J2 = 0, respectively; 
(c) and (d) are the s-wave and the d-wave solutions for J2 = 
0.05eV. (e) is the variation of the gap magnitude along the 
FS for the s-wave pairing state shown in figure (c). The inset 
in (e) is the illustration of the FS in the folded BZ (dashed 
lines). The angle 9 represents the position on the FS with 
respect to the F — M direction, a labels the inner FS (blue), 
while P labels the outer one (red). 



wavevector at which the resonance occurs is essential to 
a nonzero spin respons o 39 ' 40 . However, we note that the 
s-wave and <i-wave states are nearly degenerate around 
J2 = 0.045eV, as seen from Fig. [3] In this case, an 
s + id pairing state may be expected, which has been 
discussed in the early days of Fe pnictides 4 ^. It has been 
shown that the mixed state with s + id pairing symmetry 
does give rise to the spin resonance observed in the INS 
experiments 4 ^. On the other hand, the gap magnitude for 
the s + id pairing function is A(fe) = y/ A^(k) + A^(fc), 
where A s (fe) and Ad(k) are the gaps for the s-wave and 
the d-wave states, respectively. It is nodeless, consistent 
with the ARPES results^ii 4 .^. 



Because the spin fluctuation is dominant over the 
charge fluctuation, the pairing interaction in the spin- 
singlet channel is positive (repulsive) [see Eq. (fTTj)] and 
has a maximum around the wave vectors Q at which the 
spin fluctuation has a peak. For a repulsive pairing in- 
teraction, the most favorable SC gap must satisfy the 
condition A(fe)A(fe + Q) < in order to get the largest 
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FIG. 5. Real-space structures of the pairing gaps at different 
orbits for the s-wave and the d-wave symmetry for J2 = 0. 
The values are shown with the relative intensities of the pair- 
ing gap. d X z, d yz and d xy indicate the ti g orbits of Fe, which 
contribute mainly to the electronic states near the Fermi sur- 
face. 



eigenvalue A of the "Eliashberg" equation, as can be seen 
from Eq. (fTU|) . In the case of J2 = 0, the pairing inter- 
action comes mainly from the spin fluctuations resulting 
from the nesting of the FS. As demonstrated above in 
Fig. O (a), two sets of spin fluctuations have been ob- 
served, with peaks around q = (±n, ±0.67r), (±0.67r, ±n) 
and (±7r, ±7r) corresponding to the nesting wavevectors 
shown in FigQJb). Therefore, the favorable gaps on FS 
connected by these nesting vectors will change the signs. 
This gives rise to two possibilities, namely the s-wave 
and the d-wave symmetry, as shown in Fig. 3] (a) and 
(b). Energetically, the d-wave state here will be more 
favored because it opens a full gap around the Fermi 
pockets[Fig. Bib)]. Physically, the J2 term will induce 
the electron pairings along the NNN bonds, which gives 
rise to the pairing function of the form cos k x cosk y , i.e., 
the so-called s" 1 pairing in the case of both electron and 
hole pockets. In the absence of the hole pockets here, it 
exhibits itself as the s-wave gap function as shown in 
FigHc). 

In order to quantify the detail gap function and gain 
more understanding of the physical picture discussed 
above, we will consider the following Fourier transfer- 
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FIG. 6. Real-space structures of the pairing gaps at dif- 
ferent orbits for the s-wave and the d-wave symmetry for 
Ji — 0.05eV. Notions are the same as in Fig. [5] 



mation of the gap function: 

A^( J R) = i^e i "-«A^(g) 



(13) 



at the lowest Matsubara frequency uj n = ttT, where R 
is the real-space vector, a and /? are the orbital indices. 
We find that the inter-orbital (a ^ (3) pairing magnitudes 
are negligible in comparison with the intra-orbital (a = 
/3) pairings, suggesting that the superconducting gaps 
come mainly from the intra-orbital pairings. In addition, 
the pairing magnitudes in the e g (d z 2 and d x 2_ y 2) orbits 
are also much smaller than those in the ti g (d xz , d yz 
and d xy ) orbits, because the contribution to the FS are 
mainly from the d xz , d yz and d xy orbits. Thus, we will 
only consider the gap properties in the ti g orbits in the 
following. 

The real-space structures of the gap functions A Q( g (R) 
in the t2 g orbits up to the next nearest neighbors for 
J 2 = and J 2 = 0.05eV are plotted in Fig. E]and Fig. 
respectively. The general features are: i) the on-site par- 
ing magnitude is always zero, simply because the strong 
on-site Coulomb repulsion prohibits the on-site pairing, 
ii) The magnitudes of the d xz and the d yz orbits have only 
the Ci symmetry due to the orbital character. In the case 
of J2 = 0, the largest magnitude of the gap function ap- 
pears at the d xy orbit in the d-w&ve pairing channel. It 
is nearly twice as those at the d xz and d yz orbits in the 
same channel, and is also larger than those in the s-wave 



paring channel, so this pairing state is dominant. In this 
case, the pairing at the d xy orbit occurs only along the 
NN bonds and the pairing amplitude changes sign be- 
tween the x and y directions, so it gives rise to the usual 
d-wave gap function of the form cosfc^ — cosk y . The in- 
duction of the pairings along the NNN bonds due to the 
J2 term is obvious from Fig. [SJ But, the important fea- 
ture is that, only at the d xy orbit in the d-wave channel, 
is the NNN bond pairing not induced at all, all others 
are nearly induced equally and these pairing are larger 
than those along the NN bonds. Thus, the s-wave pair- 
ing is favored in the case of J2 = 0.05eV, which has the 
gap function of the form cos k x cos k y due to the parings 
along the NNN bonds. 

From the real-space distributions of the pairing bonds, 
we can obtain the expressions of the pairing gaps at dif- 
ferent orbits. For J2 = 0, the favorable pairing has the 
d-wave symmetry and the pairing functions are, 



A xz (k) = A [0.63cos(fc x ) - 1.26cos(/c y ) 
+ 0.14cos(A: x ) cos(k y )], 



(14) 



A yz (k) = A [1.26cos(fc a; ) - 0.63 cos(k y ) 

- 0.14cos(fc :c ) cos(fcy)], (15) 



A xy (k) = A [2. 57cos(fc E ) - 2.57cos(/c y )], 



(16) 



While, for J2 — 0.05eV, the favorable pairing has the 
s-wave symmetry and the pairing functions are, 

A xz (k) = A [0.06cos(/c 2; ) + 1.15cos(A;j / ) 

+ 2.84cos(fc x )cos(fe g )], (17) 

A yz (k) = Ao[1.15cos(fc K ) + 0.06cos(fc ?/ ) 

+ 2.84cos(fc x )cos(fc y )], (18) 

A xy {k) = A o [-0.02 cos^) - 0.02cos(fc y ) 

+ 2.90cos(A: 2; )cos(fc y )], (19) 

These pairing functions may be helpful for the experi- 
mental analysis of the pairing states. 

IV. SUMMARY 

In summary, we have investigated the spin fluctuations 
and the pairing symmetry in Aa;Fe2-j,Se2 which has only 
electron Fermi pockets, by the fluctuation exchange ap- 
proximation. Besides the on-site interactions, the next- 
nearest-neighbor antifcrromagnetic coupling J2 is also in- 
cluded. We find that both the itinerant and the localized 
natures of electrons are important to describe the recent 
experimental results of the spin fluctuations and the pair- 
ing symmetry in these materials. In particular, a small J2 
can change the pairing gap from the e?-wave symmetry to 
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the s-wave symmetry, which is consistent with the node- 
less gap observed in recent experiments. In addition, we 
have studied the real-space structures of the pairing gap 
at different orbits, from which we obtain the expressions 
for the favorable pairing gap functions. 
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